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' We address the effect of classical noise on the dynamics of quantum correlations, en- 

tanglement and quantum discord, of two non-interacting qubits initially prepared in a 
Bell state. The effect of noise is modeled by randomizing the single-qubit transition 
amplitudes. We address both static and dynamic environmental noise corresponding 
to interaction with separate common baths in either Markovian and non-Markovian 
regimes. In the Markov regime, a monotone decay of the quantum correlations is found, 
whereas for non-Markovian noise sudden death and revival phenomena may occur, de- 
pending on the characteristics of the noise. Entanglement and quantum discord show the 
same qualitative behavior for all kind of noises considered. On the other hand, we find 
that separate and common environments may play opposite roles in preserving quantum 
correlations, depending on the noise regime considered. 

Keywords: Entanglement, quantum discord, classical noise 

1. Introduction 

Quantum entanglement is a fundamental resource for quantum information pro- 



cessing, communication and exponential speed-up of some computational tasks 
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As a consequence, in the last years, there has been an increasing interest in qua n- 
tum correlat ions in various physical helds, ranging from quantum optics I 3 | 4 | 5 | 6 | ^ 
nanophysics ^™ Furthemore, recently, it was pointed out that nonclassical cor- 
relations ex ist whic h are more general, and possibly more fundamental then en- 
tanglement [Slinill] the so-called quantum discord (QD), which is defined as the 
difference between total and classical correlations in a system. Under suitable con- 
ditions QD has be en proved to be more robust than entanglement with respect to 
11 911 "31 

decoherence 1 1 I Decoherence is indeed the main threat to the possibility to fully 
exploit the potentialities of quantum correlations for the above mentioned tasks. In 
fact, the unavoidable interaction of systems with their environments induces loss 
of coherence, making the quantum parallelism, essential for quantum computation, 
ineffective. On the other hand, environment can even resume quantum correlations 
or preserve them. Thus, it is very important to analyze the effect of the various 
kinds of environmental noise on the entanglement and QD dynamics in realistic 
quantum systems, which can exhibit peculiar phenomena, such as entanglement ^31 
and QD revival. Markov noise is ascribed to environme nts with short, or rather 
instantaneous, self-correlations Non-Markovian noise is associated to envi- 
ronments with memory and may lead to the non-monotonic time dependence of 
entanglement and QD. Indeed, revival phenomena are found both for couple s of 
qubits interacting directly or indirectly in a common quantum reservoir EE1 and 
for noninteracting qubits in independent non-Markovian quantum enviro nmen ts^^. 
Non-Markovianity is also studied as a resource for quantum technology I ^U^ I. 

The effect of the environmental quantum noise on the entanglement dynamics 
between quantum systems has been interpreted in terms of the transfer of the corre- 
lations back and forth from the two-qubit system to the various parts of the global 
system. However, recent works showed that the non-monotonic time dependence 
of the amount of quantum correlations may occur in two-qubit systems under the 
local actio n of a system- unaffected environment, such as classical random external 
potentials I 20 | 21 | j) ue ^ ^ ne classical nature of the noise, in this case no back- 
action-induced correlations can be transferred from the system to the environment. 
Thus, the occurrence of entanglement revival in a non-Markovian classical envi- 
ronment is in contrast to the well-established interpretation of revivals in terms of 
system-environment quantum back action, and rises the fundamental question of 
how one could explain the effect of classical noise on quantum correlations dynamics 
in bipartite systems. 

In this work, we intend to analyze the role played by a classical noisy environ- 
ment into the entanglement and QD dynamics between two quantum particles in 
a simple system, which not only allows us to relate the intrinsic features of the 
two-particle dynamics to the revival of the correlations, but can also be of great 
relevance in various physical phenomena. Specifically, we consider two qubits non 
interacting between each other and coupled to noise in different and common envi- 
ronments. Two different kind of noise are considered: static and random telegraph 
noise (RTN). The former has also been used to describe electron transport 
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and photon propagation^ in disordered structures. The latter represents one of 
the common environmental noises affecting charge carriers in nanodevices and also 
plays a key role into the buildin g up of 1// Q noises appearing in a large number 
of solid-state systems^EIESEl] j n ^i[ s paper, we used an analytical approach to 
solve the two-qubit model, and therefore to estimate its time evolution. Numeri- 
cal techniques are also adopted to evaluate the dynamics of quantum correlations. 
Static noise is here used to simulate a non-Markovian environment. To this aim 
random time-independent terms are inserted in the off-diagonal coefficients of the 
single-qubit Hamiltonians. On the other hand, a dynamic disorder can model both 
a non-Markov environment, expressed by a slow RTN, and a Markovian noise, in 
the limit of fast RTN. In both cases single-qubits transition amplitudes are time- 
dependent and are assumed to stochastically switch between two values. 

The paper is organized as follows. The physical model adopted is described in 
Sec. 2. In Sec. 3 are presented the results for the various cases considered: static 
and dynamic noise, different and common environments, Markov and non-Markov 
regime. Conclusions and outlooks are given in Sec. 4. 

2. The Physical Model 

In this section we describe a model of two qubits, initially entangled, subject to a 
noisy classical environment. A static and a random telegraph noise are accounted for 
in different conditions, specifically local and non-local interactions between qubits 
and environments are considered. Furthermore the estimators of quantum correla- 
tions are described. 

The dynamics of the system is ruled by the Hamiltonian: 

H{t) =H A (t)®I B +lA®H B (t) , (1) 

where (t) is the single-qubit Hamiltonian defined as: 

(2) 

^A(B) and dxA(B) are t ne identity operator and the Pauli matrix of the subspace 
of the qubit A(B), respectively, s is the qubit energy in the absence of noise (here 
we assume energy degeneracy) , v is the system-environment coupling constant and 
c A(B)(t) is a random parameter related to the specific characteristics of the noise. 
This model has already be en u sed to describe a quantum walk of two not-interacting 
particles in a noisy lattice ^Sl. 

The Hamiltonian given in Eq. ([T]) is stochastic due to the randomness of the 
parameter c(t), thus leading to a stochastic time evolution of the quantum states. 
Once a choice of the noise parameter is performed, the corresponding evolution 
operator is given by U(t) — e ~ l f H (t) dt (hereafter h=l). When the latter is applied 
to the initial state the specific system dynamics is obtained. Finally, the density 
matrix describing the two qubits is evaluated by performing an average over the 
different noise configurations. The two qubits are initially prepared in the Bell state 
|$+>=1/v / 2(|0 a b ) + |1a1b)). 
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2.1. Static noise 

The static noise has already been used to study the propagation of quantum parti- 
cles i n disor dered systems, specifically in optical coupled waveguides and quantum 
walks ESl28j j n a g reemen t with these works, to model the static noise the adimen- 
sional parameters c(t) are assumed to be time-independent random variables follow- 
ing the flat probability distribution given by P(c)=l/A c for |c — co| < A c /2 and 
otherwise, cq denotes the mean value of the distribution and A c quantifies the disor- 
der of the environment. The autocorrelation function of c reads (5 c(t) 5 c(0))= A 2 ./ 12, 
and therefore its power spectrum is given by a ^-function centered on zero frequency. 
This means that the memory effects of the static noise do not vanish at any time. 
As a consequence, we can classify such a noise as non-Markovian. 

The local coupling between qubits and environment has been mimicked by as- 
suming that the noise parameters of the two qubits ca and cb are uncorrelated and 
described by two flat probability distributions each characterized by the same Cq 
and A c . On the other hand, for the case of two qubits interacting with a common 
environment we set ca=cb, which means that the same random noise parameter is 
used in both single-qubit Hamiltonians of Eq. ([2]) to model the noise effects. 

As stated above, to describe the full dynamics of the two-qubit system subject 
to the disordered environment the average over all the possible noise configurations 
is required. For the static noise, such an average is given by the integral of the time- 
evolved states each corresponding to a specific choice of the noise parameters. In 
other words, for the case of local coupling to different environments the two-qubit 
density matrix pd e (t), at time t, reads 

A c A c 

Pde(t)= / / dc A dc B P{c A )P{c B )p{cA,CB-,t) (3) 

/ A c / A c 

J C if- J c f- 

where p(cA 7 CB,t)=U(cA, cb, t)\$> + )($> + \U' f (ca, cb, t). When the two qubits are cou- 
pled to a common environment, the mixed state of the system p ce (t) is given by 

i A c 

p ce (t) = / dcP(c)p(c,t) (4) 

/ A c 
J CO f- 

where p(c,t)=U(c,t)\$ + )($ + \W (c,t). The explicit form of pdeif) and p ce (t) is pre- 
sented in Sec. [3] 

2.2. Random telegraph noise 

The other kind of noise we examine is the RTN, which is able to de scribe a number 
of typical phenomena affecting solid-state devices on the nanoscale l ^4 | 25 | jj ere c (£^ 
is assumed to flip randomly between the values -1 and 1 at rate 7. The autocorre- 
lation function falls off exponentially as (8c(t)8c(Q))—e^ 2lt and its power spectrum 
exhibits the well-known Lorcntzian shape 47/ (uJ 2 + 47 2 ). Two different regimes of 
the decay of quantum correlations arc identified, according to the ratio between the 
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system-environment coupling v and the switching rate 7 of RTN^^. For 1//7 -C 1 
the Markovian regime is found, while for » 1 the non-Markovian behavior is 
obtained. 

As discussed by Abel et a/.^l the dichotomic stochastic behavior of c(t) induces 
a random phase factor ipA(B)(t) = ~ v J dt'cA(B)(t') in the evolution of the single- 
qubit states. The two-qubit density matrix is given by the average over the random 
phase factor, namely 

p de {t) = ({p{<p A {t),<p B {t))) VA ) ^ (5) 

for local coupling qubit-environment, and 

Pce (t) = { P {^{t))) v (6) 

for coupling of the two qubits with a common environment. Also in this case, the 
explicit forms of pde{t) and p ce {t) are given in Sec.[3l 



2.3. Estimators of quantum correlations 

Here, we briefly review the estimators used to quantify entanglement and quantum 
discord in our system. 

Entanglement is evaluated in terms of negativity SOI defined as: 



(7) 



where A~~ are the negative eigenvalues of the partial transpose of the density matrix 
of the bipartite system. Negativity ranges from zero, for separable states, to one, 
for maximally entangled states. 

QD is the difference between the total and the classical correlations in a system 
so representing its degree of quantumness ™ It is defined as 

Q = I-C 1 (8) 

where I is the quantum mutual information which gives a measure of the total 
amount of correlations of the system. It is defined as: 

l = S(p A ) + S(p B )-S(p), (9) 

where S is the von Neumann entropy and p A ^ BS> indicates the reduced density ma- 
trix of the subsystem A(B). C denotes the measurement-induced quantum mutual 
information, namely the classical correlations. The latter reads 

C = m a *{S(p A )-S(p\{B k })}, (10) 

{B k \ 

with S(p\{Bi c }) indicating the quantum conditional entropy with respect to the set 
of measurements {Bk} performed locally on subsystem B. Usually to compute QD 
is not an easy task, since it involves a maximization procedure. But i n th e case of 
two-qubit X states an analytical expression for Q was derived by Luo EH. 
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3. Results 

In this Section, we present the explicit forms, as a function of time, of the various 
mixed states of the two qubits, and the time behavior of entanglement and QD 
for all the considered physical configurations. Explicit derivation for the average 
density matrices is shown in | Appendix A| 



3.1. Static noise 



When the two qubits are affected by static noise, their density matrices, obtained 
from Eqs. ([3]) and (|4j) , take the form: 



Pde(ce) (t) 



/ 4 + a de(ce) —Pde(ce) ~Pde(ce) \ + a de(ce) \ 
Pde(ce) 4 ^de(ce) 4 ^de(ce) ftde(ce) 



^de(ce) 4 
Pde(ce) 4 ^de(ce) 4 &de(c 



(11) 



V \ + a de(ce) —fide(ce) —^de(ce) \ + &de(ce) I 



where 



1 



ade 



(2A c ut) 2 
1 



cos (4cqi4) sin 2 (A c i>t) and (3, 



de 



8A c vt 



cos (icQft) sin (2A c vt) and (3 C 



(2A C ^) : 



8A c vt 



; sin (4co^£) sin 2 [A c vt) 
sin (AcQvt) sin {2A c vt). 



It should be noticed that, for both local and non-local qubit-environment coupling, 
at sufficiently long times, the matrices take an X form where all non-vanishing 
coefficients are equal. This means that the asymptotic two-qubit state is given by 
a statistical mixture of two Bell states, specifically 1/2 (|$+)($+| + 

Performing the calculation starting from the definition given in Eq. ([7]), the 
negativity can be be expressed as 



Mde(t) = 



sin (A c vt) 
A r vt 



and JVce{t) = 



sin {2A c vt) 
2A r M 



(12) 



Such a result clearly shows a non-monotonic time decay of entanglement as expected 
from the non-Markovian nature of the noise. Peculiar phenomena as sudden death 
and revival are observed (see left panel of Fig. [1]). Furthermore, we find that in 
the case of common environment entanglement is more robust with respect to the 
case of different environments. Indeed, when the two qubits are coupled with a 
common environment, the latter can be interpreted as a sort of interaction mediator 
between the qubits themselves. Such an interaction somehow contributes to build 
up quantum correlations even if decohering effects of the environment are still 
dominant and lead to a power-like decaying profile of entanglement. An analogous 
behavior is shown by QD, as displayed in the right panel of Fig. [TJ Note that, 
unlike negativity, QD has been evaluated numerically by means of a maximization 
algorithm. 
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2 4 6 8 10 12 14 16 18 2 4 6 8 10 

time (arb. unit) time (arb. unit) 



Fig. 1. Time evolution of negativity (left panel) and QD (right panel) for two qubits subject to a 
static noise for different and common environments. To highlight the revivals of QD, its behavior 
is reported also in logarithmic scale in the inset of the right panel. 



3.2. Random telegraph noise 

For the case of a RTN affecting the dynamics of two qubits, the time evolution of 
the mixed state is evaluated from the average over the random phase factor <p(i). 
After a straightforward calculation, we obtain 



1 + A 



de(c 



e(ce){t) - J 



\1 + A 



de(c 









^de(ce) 


1 - Ade(ce 


Arfe(ce) 


1 - A de ( ce 









1 



Ade(ce) \ 









1 + A 



(13) 



de(ce 



J 



with Kd e =D 2 v (t) and A ce =£>4„(£). As shown elsewhere the function D(t) cor- 
responds to the average phase factor (e iv ^>) and can be expressed as 



,- 7 t 



D mv (t) 



-it 



cosh (8 mv t) + sinh (S mv t) 
cos (8 mi ,t) + -f 1 - sin (S mv t) 



for 
for 



7 > mv 
7 < mv 



(14) 



where S mu =^y\j 2 — (mv) 2 \ with m £ {2,4}. Unlike the static noise, the two-qubit 
density matrix assumes an X form both for different and common environment at 
any time t. It is worth noting that, for both RTN configurations, the long time 
asymptotic form of the states obtained is the same found for the static noise. Using 
again the definition of Eq. (J7J, negativity reads: 



N de (t) = D\ v (t) and Af ce {t) = \D iv (t)\. 



(15) 



In this case, due to the X form of the states, it is possible to give the analytical 

lOI I 

expressions for QD by following the Luo approach . They read: 

Qde(t) = \ [(1 + D\ v ) log 2 (1 + D\ v ) + (1 - D\ v ) log 2 (1 - D\ v )\ 



Q ce (t) = - [(1 + D 4v ) log 2 (1 + D 4u ) + (1 - D iv ) log 2 (1 - D iu )} . (16) 
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The dynamics of Af and Q is shown in Fig. [2] In agreement with the outcomes 

IOOI 

of previous works'^ 21 , in the Markovian regime both entanglement and QD exhibits 
an exponential decay with time. On the other hand, in the non-Markovian regime 
the amount of quantum correlation are a damped oscillating function of time, thus 
displaying sudden death and revival phenomena. Finally, it should be noticed that 
our calculations show that, in the Markovian regime, the quantum correlations for 
the case of coupling of two qubits with the common environment are weaker with 
respect to the case of local qubit-environment interaction. Such a behavior is in 
contrast with the one found for non-Markovian noises and can be ascribed to the 
fact that indirect qubit-qubit interaction is here more effective in destroying both 
entanglement and QD. 




Fig. 2. Top panels: Time evolution of negativity (left) and QD (right) for two qubits subject to a 
RTN for different and common environments when 1//7 = 0.2 (Markovian regime). Bottom panels: 
Time evolution of negativity (left) and QD (right) for two qubits subject to a RTN for different 
and common environments when v/j = 5 (non-Markovian regime). To highlight the revivals of 
QD, its behavior is reported also in logarithmic scale in the inset of the right panel. 



4. Conclusions 

In this paper we have analyzed the effects of both independent and common noisy 
environments on the dynamics of quantum correlations of two non-interacting 
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qubits. Negativity and QD have been computed, using both analytical and numer- 
ical techniques. Static and slow RTN are classified as non-Markovian noises, while 
the fast RTN mimics a Markovian environment. We showed that, starting from 
a maximally entangled state, the quantum correlations display different decaying 
behaviors, depending on the nature of the considered noise. In particular, Marko- 
vian environments lead to a monotonic decay, while in non-Markov regimes sudden 
death and r evival phenomena are present, in agreement with previous results in 

loniooi 

the literatur e 1 '. Specifically, here we ascribe the decay of revivals obtained in 
the static noise to the continuous nature of the random noise parameter. In fact 
our analysis clearly shows the suppression of quantum correlations at long times 
unlike previous investigations where an environmental noise mimicked by means of 
dichotomic time-independent noise parameter leads to non decaying entanglement 
revivalJ^M 

For both static noise and slow RTN, our results highlight that a common environ- 
ment preserves better quantum correlations. It is worth noting that the opposite is 
found for qubits subject to a fast RTN, where the effect of a common noise results in 
a faster decay of correlations with respect to the case of independent environments. 

An interesting generalization of our model is the possibility to sum up many 
RTNs, each with a specific switching rate, in order to obtain 1/ f a noise. This noise 

IOOI 

is ubiquitous in solid state devices^ and can be simulated considering a linear 
superposition of different RTNs, with an appropriate distribution of their switching 
rates. Furthermore, our study of the effects of static and dynamic noise on the 
evolution of quantum correlations in two-qubit models, has a natural and quite 
promising generalization in the case of systems with a larger number of degrees of 
freedom. For instance the model here studied can be used to extend, to t he p resence 
of noise, previous works on quantum walks on a one-dimensional lattice^. In this 
case the environmental noise would be a consequence of the lattice disorder. Work 
along these lines is in progress and results will be reported elsewhere. 
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Appendix A. Explicit derivation of the two-qubit density matrix 

Given the Hamiltonian in Eqs. ([T|) and ©, the associated evolution operator reads: 



Ude(ce) =e 



-2i0t 




(A.l) 
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where 



A de = cos(tp A (t) ) cos((p B (t)) 
B de = cos(<pA(t))sin(<p B (t)) 
C de = sm(ip A (t)) cos(ip B (t)) 
D de = sm((p A (t))sm((p B (t)) 



A ce = cos%(i)) 

B ce = cos(ip(t)) s'm((p(t)) 

C ce = B ce = cos((p(t)) s'm(ip(t)) 

D ce = sm 2 (ip(t)), 



with the phase ip(t) — ct in the case of static noise and ip(t) = —v J* c(t')dt' for 
the RTN. Applying the evolution operator to the initial state po — |$ + )($ + |, it is 
possible to obtain the time-evolved density matrix of the system corresponding to 
a specific choice of the random parameter: 



e(ce) 



/ -^de{ce) -^de(ce) ^de(ce) ^de(ce) \ 



de(ce) ^de(ce) 



Cdi 



\ A 



de(c 



(ce) 

Cde(ce) Ccie(ce) 



de(ce) -^de(ce) ^de{ce) 



B d 



B 

-B 
A d( 



de(ce) 
de(ce) 
(ce) 



(A.2) 



where: 



A-de(ce) — (^-rie(ce) — ^de(ce)) 2 

Bde(ce) — i{B de (ce) + Cde(ce) ) ( Aie(ce) — Z?d e (ce)) 
Cd 



(A.3) 



^de(ce) — (B de (ce) +C(2e(ce)) 

The average of the expression IA.2I over the possible values of the noise parameters 
gives the state of the system at time t and corresponds to the integral 



d(p(. . . )p(cp, t) 



(A.4) 



where for the static noise the phase distribution p(w, t) — P(c)/t where P(c) is the 
flat distribution given in Sec. 12.11 while for the RTN, it takes the forr 



p((p,t) = -e" 7 * ([% + vt) + <% - vt)] + -[e(cp + vt)+Q(<p- vt)]) 
2 1 v > 

=-L + I (jty/l - (tp/vty 



(A.5) 



Ik(x) is the modihed Bessel function and 9(x) is the Heaviside step function. After 
performing the integration, the average density matrices take the form of Eqs. (jlip 
and (fT3|) for the static noise and the RTN respectively. 
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